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Abstract 

The exponential scaling of isotopic yields is investigated for sources of different sizes over a broad 
range of excitation energies and freeze-out volumes, in both primary and asymptotic stages of the 
decay in the framework of a microcanonical multifragmentation model. It was found that the 
scaling parameters have a strong dependence on the considered pair of equilibrated sources and 
excitation energy and are affected by the secondary particle emission of the break-up fragments. No 
significant influence of the freeze-out volume on the considered isotopic ratios has been observed. 
Deviations of microcanonical results from grandcanonical expectations are discussed. 

PACS numbers: 25.70.PqMultifragment emission and correlations, 24.10.PaThermal and statistical models 
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I. INTRODUCTION 



Many nuclear physics experiments realized in the last years revealed an interesting scaling 
property of the isotopic yield ratios obtained from the disintegration of equilibrated sources 
of similar sizes, with close values of excitation energy per nucleon and temperature but 
different isospin values as a function of the isotopic composition of the emitted cluster 



This property has been called isoscaling and can be expressed mathematically by the 



formula, 



R 21 (Z,N) = Y 2 (Z,N)/Y 1 (Z,N) = Cexp(aN + pZ), (1) 



where Yi(Z, N) denotes the yield of the isotope (Z, N) obtained from the decay of the excited 
nuclear system "z" and C, a and (3 are scaling parameters. 

If one classifies the nuclear reactions which manifest isoscaling function of the projectile 
energy, one can say that isoscaling has been observed in reactions induced by projectiles 
whose energies range from a few MeV/nucleon to several GeV/nucleon. In terms of reaction 
mechanisms, isoscaling has been obtained in evaporation reactions (e.g. 4 He + 116 Sn, 124 Sn 
at 50 MeV/nucleon bombarding energy deep inelastic reactions (e.g. 16 + 232 Th, 16 
+ 197 Au at 8.6 MeV/nucleon 4], 86 Kr + 112 ' 124 Sn and 86 Kr + 58 < 64 Ni at 25 MeV/nucleon 
[^j]) and, recently, in fission reactions (e.g. n + 233 > 238 U at 14 MeV/nucleon jfjj). The 
biggest amount of experimental isoscaling data comes from multifragmentation reactions 



and the great interest in these studies was motivated by the possibility 
to infer information on the nuclear equation of state, and in particular on the asymmetry 
term, from the isoscaling parameters [a, . 

It is interesting to specify at this point that, relatively easy to understand for nuclear 
reactions involving equilibrated compound nuclear systems isoscaling is compatible also 
with dynamical fragment formation as was recently demonstrated in Ref. [12I 

The aim of the present paper is to investigate isoscaling in multifragmentation reactions 
and more precisely the results produced within a microcanonical framework. In order to learn 
how sensitive the isoscaling parameters are with respect to the observables characterizing 
the source's state, we considered sources whose nucleon number varies from 40 to 200, with 
excitation energies between 2 and 15 MeV/nucleon, freeze-out volumes ranging from 3Vq 
to IOVq in both primary and asymptotic stages of the decay. The relevance of such study 
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consists in the fact that the most suitable statistical ensemble to describe the decay of 
an isolated excited nuclear system counting at most several hundreds of nucleons is the 
microcanonical ensemble. The importance of using models which can mimic as well as 
possible the real physical case is demonstrated by the fact that the best description of the 
experimental data was obtained by the predictions of microcanonical multifragmentation 



models 
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Section |H] makes a review of the grandcanonical argumentation of isoscaling in nuclear 
multifragmentation. Section IIHI presents microcanonical results on isoscaling obtained for 
different mass systems with excitation energies ranging from 2 to 15 MeV/nucleon at differ- 
ent freeze-out volumes in both primary and asymptotic stages of the decay. The sensitivity 
of the isoscaling parameters to observables characterizing the state of the source is dis- 
cussed in detail. Deviations of microcanonical results from grandcanonical expectations are 
commented. The conclusions are drawn in Section HVl 



II. ISOSCALING IN NUCLEAR MULTIFRAGMENTATION 

Isoscaling in multifragmentation was for the first time evidenced for the light emitted frag- 
ments (A < 18) obtained in the reactions 112 Sn+ 112 Sn and 124 Sn+ 124 Sn at 50 MeV/nucleon 
bombarding energy [l| and confirmed by many other reactions induced by relativistic light 
oarticles (e.g. p, d, a+ 112 Sn and 124 Sn at incident energies between 660 MeV and 15.3 GeV 
^) or heavy ions ( 124 Sn+ 64 Ni and 112 Sn+ 58 Ni at 35 MeV/nucleon 0). 

The theoretical interpretation of isoscaling in the multifragmentation energy domain was 
given using a grandcanonical ensemble |7[. It was demonstrated that in the grandcanon- 
ical limit the yield of a primary fragment (Z, N) obtained from the decay of a thermally 
equilibrated system characterized by the temperature T can be written as 

Y(Z, N) oc exp(B(Z, N)/T) exp(iV> n /T + Zfi p /T), (2) 

where B(Z,N) is the binding energy of the emitted fragment (Z,N) and \x n and \i p are 
the neutron and proton chemical potentials. Taking into account that the two equilibrated 
systems have the same temperature, excitation energy per nucleon and similar sizes the 
proportionality factors from Eq. (j2J have been considered practically identical in the two 
cases. Thus, one obtains for the ratio of the yields of any isotope resulting from the decay 
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of the two sources the expression of Eq. (JIJ) where, 



a = A/i n /T, (3) 

and 

P = A/Up/T, (4) 

where Afi n{p) = fi 2n ( P ) - Hin(p)- 

Despite the fact that the suitability of a grandcanonical approximation is hard to be 
justified for systems having at most few hundreds particles, the success of the above formulas 
was due to their simplicity and to the fact that no contradiction came from experimental 
data or various statistical multifragmentation models (e.g. microcanonical Q| and canonical 
[lsj ] versions of SMM, Expanding Emitting Source ^) currently used in data analyzes 

Concerning the slope parameters a and /3, several values have been reported in literature. 
This fact is expected because is known that particle formation probabilities have explicit 
dependence of all observables which characterize an excited source (mass A, charge Z, 
excitation energy E ex , freeze-out volume V). It is also known that the sequential particle 
emission from the excited primary fragments modifies the isotopic yields. This evaporation 
process may also affect the values of a and /3. 



III. ISOSCALING FROM A MICROCANONICAL PERSPECTIVE 

Despite the fact that, contrary to the grandcanonical case, no microcanonical model pro- 
vides analytical expressions for cluster yields, the occurrence of isoscaling in a microcanonical 
model can be understood starting from the shape of the isotopic yields corresponding to a 
given element Thus, it has been observed that charge distributions of fragments with 
fixed values of the neutron number (Y(Z, N)\n) and N distributions of fragments with 
fixed values of the proton number (Y(Z, N)\z) can be approximated by Gaussian functions. 
The occurrence of such Gaussian-like distributions within a microcanonical model can be 
regarded as a natural consequence of complete equilibration. Using the generic notation 
Y(x, y)\ y for the distribution of the observable x when y is kept fixed (where x = N, Z and 
y = Z,N), one can write, 



Y(x,y)\ y = Cexp(-(x - x me6 {y)f /2a{y) 2 ). (5) 
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The mean value x me d (y) gives the isotope produced with the largest probability, while the 
variance a(y) is a measure of how uniform is the population of fragments having different 
isospin values. The statistical meaning of these quantities involves an obvious dependence 
on observables determining the ensemble's state. Since we deal with a microcanonical ap- 
proximation, these observables are the mass, volume and excitation energy of the multifrag- 
menting source. 

Using Eq. the ratio of isotopic yields writes, 



D / . N Y 2(x,y)\ t 

R2i{x, y = const. 



Yi(x,y)\ 



exp 



x 2 ( 1 \ ( x med2 (y) x medl (y)\ ( x med2 (y) 2 x medl (y) 2 



(6) 



2 \a 2 (yf - aM 2 J Wl/) 2 °i(y) 2 J V ^(y) 2 2a x {yf 
Taking into account that the two considered sources have similar sizes, temperatures and 
freeze-out volumes, the isotopic distributions corresponding to any cluster are expected to 
have equal variances (<Ti(y) = cr 2 (y) = cr(y)). Thus, in Eq. © the quadratic term in x 
vanishes leading to a linear dependence of the yield ratio on the observable x of the emitted 
cluster. 

From Eq. © one can express the a and f3 slope parameters in terms of mean values and 
variances of isotopic distributions, 

a{Z) = -i- (N med2 (Z) - N medl (Z)) , (7) 
a(Z) 2 

m = (^WiV) - Z medl (iV)) • (8) 

From Eqs. (jZJ) and (jHJ) results that a and j3 may depend on the size of the emitted 
cluster. A necessary condition to obtain slope parameters independent of the size of the 
emitted cluster as in Eq. Q is to have a constant value for the ratio of the shift between 
the isotopic distributions produced by the two sources and the square of the distributions' 
variance. As we shall see in the following sections this hold with good approximation for the 
light clusters emitted by relatively large equilibrated sources but does not hold for heavy 
fragments emitted by heavy sources and clusters emitted by light multifragmenting sources. 

In the following we shall investigate isoscaling over a large range of masses, excitation 
energies and freeze-out volumes of the multifragmenting sources in order to determine the 



dependence of the slope parameters a and 3 on the source parameters, the influence of 
the secondary decays on a and (3 and, in the case of small sources, the dependence of a 
and 3 on the size of the emitted cluster. To serve this goal we use the micro canonical 
■nu.ti^entation.ode.p.esentedindetai.inR^.Efi. 

A. Model overview 

The micro canonical multifragmentation model Q,!^ has two distinct stages: the break- 
up stage and the sequential particle evaporation stage. 

The most important part of the model is the so called break-up stage which aims to 
describe the explosion of the equilibrated nuclear source. In order to mimic as well as possible 
the physical situation of a small isolated system with fixed excitation energy, the model tries 
to obey sharply to microcanonical constrains. Thus, the fixed observables characterizing 
the state of the equilibrated source are mass, charge, total energy E, total momentum P 
(=0 in the center of mass (cm.) frame) and total angular momentum L (=0 for non- 
rotating systems). The freeze-out volume V is determined by the spherical container in 
which fragments are generated. 

The standard version of the model allows treatment of freeze-out volume according to 



two different scenarios [16]: hard spherical nonoverlapping fragments 
container and spherical container with free volume parameterization 



P 



aced in a spherical 



211 ] . For a faster sim- 



ulation we considered the second case. The mathematical expression of the free volume is: 



N 



Vfr~ = Y[(V-i-±), (9) 

i=l J 

where Nf r is the number of fragments in a given event and Vq is the volume of the nuclear 
system at normal density. 

The price to pay for a model aiming to respect all microcanonical constrains is the 
impossibility to obtain formulas analytically tractable. The key quantity of the model is the 
statistical weight of a configuration C, 



w ° « 7v^ g j x r( 3/2(7v c - 2)) {mA) w ■ (10) 



6 



where Nc is the number of fragments corresponding to configuration C, Q is the accesible 
volume, I is inertial tensor of the system and K is the available kinetic energy. The index 
n denotes the fragments in each configuration. 

The average value of any observable Q is calculated via a Metropolis type simulation 
according to the expression, 

Because the observables studied in the present paper are functions of isotopic yields they 
are strongly dependent on all specific details of the model (binding energy parameterization, 
allowed isospin asymmetry of primary fragments, employed level density formula, etc.). The 
)reak-up fragments allowed to be formed are all isotopes included in the mass table of Ref. 



221- 

The fragments' binding energy was calculated according to the liquid-drop parameteri- 
zation: 



B(A,Z) = 15.4941(1 - 1.7826/ 2 )A 

- 17.9439(1 - 1.7826/ 2 )A 2/3 

- 0.7053Z 2 A~ 1/3 + 1.1530Z 2 /A, (12) 

where I — (A — 2Z) /A. 

Fragments with mass smaller than 4 are considered without internal excitation energy; 
larger fragments are allowed to carry an excitation energy (e) upper limited by the binding 
energy according to the following level density formula: 



= 12q i/4.5/4 exp(2v^) exp(-e/r), (13) 
with a = 0.1UA + 0.098A 2 / 3 MeV" 1 Q and r=9 MeV. The factor exp(-e/r) takes into 
account the decrease of the excited level lifetime with the increasing excitation energy. 

The excited primary fragments obtained in the break-up stage are allowed to decay by 
sequential particle emission (second stage of the model) as in Ref. flfil ]. The range of the 
evaporated particles is considered up to A = 16. As already stated in literature j^J, |3] by 
modifying isotopic yields, secondary particle production mechanisms may alter isoscaling 
parameters. 



Despite the fact that the employed parameterizations are expected to affect the specific 
values of a and (3 they do not affect the general behavior of isoscaling parameters as a 
function of observables characterizing the sources' state such as the conclusions of the study 
are general and valid. 

B. General microcanonical results 

In order to verify to what extend isoscaling holds within a microcanonical framework 
and, if this is the case, to make a systematic study of all source parameters which might 
have an influence on the isoscaling parameters we considered pairs of equilibrated sources 
with masses ranging from A = 40 to A = 200, excitation energies ranging from 2 to 15 
MeV/nucleon, freeze-out volumes ranging from 3Vo to 10Vo in both primary and asymptotic 
stages of the reaction. 

For all considered cases we obtained isoscaling in the sense that log R21 (Z, N) shows an 
almost perfect linear dependence as a function of N and Z. 

As general observations, one may say that the quality of the isotopic scaling is better in 
the break-up stage than in the asymptotic stage of the decay and the maximum size of the 
emitted cluster for which isoscaling still holds decreases with increasing excitation energy. 
We think that the more modest quality of the isoscaling in the asymptotic stage is an artifact 
of the simplified manner in which the secondary decays have been implemented. We adopt 
the widely used convention to denote with the index "2" the more neutron rich system and 
with the index " 1" the more neutron poor system. In this situation the value of a is always 
positive because more neutron rich clusters will be produced by the neutron richer source 
and the value of (3 is always negative. To avoid possible effects of different magnitudes of 
Coulomb interaction on isotopic distributions, the considered pairs of sources have the same 
proton number and different mass numbers. 

To illustrate how isoscaling looks like in the case of this microcanonical model, in Fig. 
1 we present the isotopic yield ratios as a function of N and Z for the pair of equilibrated 
sources (200, 80) and (185, 80) at 5 MeV/nucleon excitation energy and a freeze-out volume 6 
times larger than the nuclear volume at normal density. As results from the figure, isoscaling 
holds in the strict sense of Eq. Q namely linear behavior of log R21 as a function of N and 
Z and constant value of C. Special attention should be drawn to the fact that in Fig. Q 
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FIG. 1: Isotopic yield ratios are plotted in logarithmic scale as a function of N (upper panels) and Z 
(lower panels) for two equilibrated sources (200, 80) and (185, 80) with the excitation energy E ex =5 
MeV/nucleon and the freeze-out volume V = 6Vb in both primary (left panels) and asymptotic 
stages of the decay (right panels). The considered sequence of N(Z) is 1,2,3,4,5,6,8,10,12,14,18. 

distances from the points corresponding to different values of Z(N) is not equal because the 
considered sequence of Z(N) does not contain exclusively consecutive numbers (see legend). 

It is interesting to notice that a linear behavior of logit^i is obtained also for emitted 
clusters much heavier than those for which the experimental analysis has been possible. 
As a technical detail, we stress that despite for a given element the linear behavior of 
log #21 (-^0 holds for almost all obtained isotopes, for the fitting procedure only the most 
stable fragments have been selected. 

C. Dependence of a and (3 on excitation energy 

Since, as known from the early days of multifragmentation, the excitation energy of the 
nuclear source induces strong modifications on all fragment multiplicities, it is expected to 
affect also isotopic yield ratios. 
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In Fig. Elare represented a and (3 for the equilibrated sources (200, 80) and (185, 80), with 
V = QVq as a function of the excitation energy. It must be specified that, if not mentioned 
in a different manner, along this paper a and (3 are calculated such as to obtain the best fit 
for fragments with 2 < Z < 9 and respectively 2 < iV < 9. This choice is made in order to 
realize an analysis as much as possible similar to the way in which the slope parameters have 
been obtained from experimental data and to have the same amount of considered data for 
all excitation energies and all considered source sizes. 

One can see that both parameters have a monotonic dependence on the excitation energy, 
their absolute value decreasing with the temperature. Thus, increasing the excitation energy 
from 2 MeV/nucleon to 15 MeV/nucleon a diminishes by a factor of 2.5 and \/3\ diminishes by 
a factor of 3.5 This effect can be understood having in mind that an increase in excitation 
energy will result in a more uniform population of the isotopes corresponding to a given 
element thus washing out the effect of the isotopic composition of the source. In terms of Eqs. 
((Zj) and (jEJ), an increasing excitation energy will lead to a shift of the Y(x,y)\ y distributions 
toward smaller values of x me d such as (x me d2 — ^medi) slightly diminishes together with an 
enhancement of the distributions' variance a. Since both numerator and denominator from 
Eqs. (J7J) and (jHJ) act in the same sense, the decrease of a and with the excitation energy 
of the sources is obvious. 

These results are in qualitative agreement with those of Ref. B llo|. 



D. Effect of secondary decays on a and (3 slope parameters 

An important problem that manifests whenever we need to infer the break-up information 
from the available experimental multifragmentation data is the reconstruction of primary 
observables affected by the secondary particle emission. This problem may be even more 
important when one operates with very sensitive quantities like the slope parameters a 
and (3. It was stated in the literature that when the two considered sources have almost 
identical parameters excepting the isospin, the effects of sequential evaporation cancel one 
another leading to small modifications of the break-up stage results Q . On the other hand, 
estimations performed with accurate models like SMM ;14] indicate that modifications of 
the slope parameters due to evaporation may manifest with different magnitudes depending 
on how sequential evaporation is implemented j^, 11 ] . 
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FIG. 2: Slope parameters a and [3 as a function of excitation energy for primary (open symbols) 
and asymptotic (full symbols) stages of the decay. The considered equilibrated sources are (200, 
80) and (185, 80). The freeze-out volume is V = 6Vq- 

To check the evaporation effect on isoscaling within the present model in Fig. El are 
plotted a and (3 in the break-up and asymptotic stages of the decay for the equilibrated 
sources (200, 80) and (185, 80), with V = 6Vo as a function of the excitation energy. As 
displayed by the figure, over the entire energy range the absolute values of a and (3 are larger 
in the asymptotic stage of the decay in comparison to the break-up stage. 

To get an insight on what happens during the sequential particle emission from primary 
excited fragments is necessary to focus on the evaporation process. Given the fact that for all 
excitation energies neutrons are emitted with the most important probability, the final cold 
fragments are more symmetric than their break-up ancestors. This means that the isotopic 
distribution of a given Z shifts toward smaller values of N me d. Since the amplitude of this 
shift is larger for the neutron richer source, it means that the value (iV me d2 {Z) — N me( n(Z)) 
in the asymptotic stage is smaller than the one corresponding to the break-up stage. This 
apparent trend of isoscaling parameters to decrease after particle evaporation is annihilated 
by a strong narrowing of isotopic distributions (the more asymmetric nuclei have small 
survival probabilities after secondary decays) such as in the end the values of a and \/3\ 
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FIG. 3: Effect of secondary decays on the slope parameters a and (3 function of excitation energy. 
The considered equilibrated sources are (200, 80) and (185, 80) with V = 2>Vq. 

increase with respect to their break-up values. This result is counterintuitive because at a 
first glance one could expect secondary decays to wash out the isotopic differences of the 
sources. 

From Fig. E] is obvious that secondary decays affect more a than (3. This result can 
be understood taking into account that the decrease of the width of Y(Z, N)\z distribu- 
tions caused by neutron evaporation is more important than the narrowing of Y(Z, N)\n 
distributions caused by less important charged particle emission. 

One can see from Fig. El that the difference between the asymptotic and break-up values 
increases up to 5 MeV/ nucleon excitation energy and then stays roughly constant, but taking 
into account that the increase in excitation energy results in smaller slope parameters, this 
means that the relative influence of the secondary decays on a and (3 grows with energy. To 
illustrate the relative modifications brought by the secondary decays to the slope parameters 
in Fig. El is plotted the deviation between asymptotic and primary values relative to the 
break-up results versus excitation energy. The increasing rate is roughly constant for whole 
considered domain. 



12 



E. Slope dependence on the freeze-out volume 



To identify a possible freeze-out volume effect on the slope parameters, isoscaling was 
analyzed for the previously mentioned pair of sources ((200, 80) and (185, 80)) at excitation 
energies increasing from 2 to 15 MeV/nucleon at V = 3V , V = 6V , V = 10V in both 
primary and asymptotic stages of the decay. Despite the fact that all isotopic distributions 
present a strong dependence on the volume of the excited source, for the considered cases 
no significant freeze-out volume dependence of a and (3 was identified. 

The impossibility to get an analytic expression of isotopic yields within a microcanonical 
model makes very hard to estimate the effect of freeze-out volume (or any other observable 
characterizing the state of the source) on isotopic ratios. However, if true, the independence 
of the slope parameters on the freeze-out volume is of particular importance for experiments 
aiming to extract the symmetry term of the binding energy from isotopic yields. Indeed, 
no matter is the method for selecting collisions, is hard to believe that the ensemble of 
equilibrated nuclear systems formed in heavy ion reactions correspond to one value of the 
freeze-out volume, but rather to a distribution. Since there is no method to select fragmen- 
tation events according to their volume, volume dependent a and (3 would make impossible 
the extraction of the asymmetry term. 



Ref. 



Ill ] makes use of a canonical version of SMM ^| and reaches the conclusion that 



for all considered temperatures (T=4, 5, 6 MeV) a increases monotonically with the freeze- 
out density: by increasing p/po from 0.1 to 1 a increases by a factor of two. The opposite 



conclusions of Ref. 11] and the present study are produced by the in principle different basic 
hypothesis of each model. Taking into account the implications of the volume dependence 
of isoscaling parameters, a definite answer to this problem would be of much interest. 



F. Slope dependence on the considered sources 

All the above discussion suggests a strong dependence of the fragment production on the 
parameters characterizing the equilibrated state of the sources. A natural question to rise 
at this point is how sensitive are the slope parameters on the considered sources once one 
keeps fixed the atomic number, the freeze-out volume and the excitation energy. 

To answer this question a and (3 have been calculated for different pairs of sources (chosen 
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FIG. 4: Slope parameters a and (3 as a function of the difference between the binding energies of 
the equilibrated sources. Open symbols correspond to break-up isotopic yields and full symbols 
correspond to isotopic yields after sequential particle emission. The considered equilibrated sources 
are characterized by Z = 80, E ex = 5 MeV/nucleon and V = 3Vq. 

from the set (180,80), (185,80), (190,80) and (200,80)) at E ex = 5 MeV/nucleon and V = 
3Vo. The result we obtained is that the slope parameters are not dictated by the difference 
of isospin asymmetry between the two sources, but by the difference between the binding 
energy of the sources. To illustrate this result we plot in Fig. HJthe monotonic dependence 
of a and (3 as a function of the difference between the binding energies of the equilibrated 
sources. 



This result is in qualitative agreement with some experimental data. Ref. |lj| reports 
a modification of the absolute values of a and (3 by a factor of two when calculated using 
( 112 Sn+ 124 Sn/ 112 Sn+ 112 Sn) nuclear reactions instead of ( 124 Sn+ 124 Sn/ 112 Sn+ 112 Sn) at 50 
MeV/nucleon. Similar results were obtained recently studying isoscaling in ( 58 Fe+ 58 Ni/ 
58 Ni+ 58 Ni) and ( 58 Fe+ 58 Fe/ 58 Ni+ 58 Ni) reactions at 30, 40 and 47 MeV/nucleon Q- 

Contrary to the cases of a grandcanonical model or a sequential particle emission model 
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where one can easily express isotopic multiplicities as a function of the binding energy 
of the parent nucleus, within a microcanonical model the effect of binding energy can be 
estimated only in an intuitive way. Thus, more similar are the two considered equilibrated 
sources from the point of view of their binding energies, closer are the values of the total 
available energy of each system for a given value of the excitation energy. This will lead to 
similar populations of the configuration space of each system. In this case, all observables 
(including isotopic yields) are expected to have close values. Close values of isotopic yields 
will obviously produce small values of a and (3. 



G. Slope dependence on the source size 

The only parameter characterizing the equilibrated source which was kept fixed thorough 
the above discussion is the source's size. To check the size effect of isoscaling parameters 
we choose for fragmenting sources two pairs of smaller nuclei, namely ((112, 50), (124, 50)) 
and ((112, 50), (119, 50)) with E ex = 5 MeV/nucleon and V = 3V . The values of the slope 
parameters obtained in this case (((112, 50), (124, 50)): 0^=0.37, a as =0.50, ^=-0.49 and 
/3 as =-0.55; ((112, 50), (119, 50)): a fefc =0.23, a as =0.32, (3 bk =-0.30 and /3 as =-0.34) are very 
different from the values obtained for the Z = 80 sources and suggest a strong dependence 
of the isoscaling parameters on source size. 

This result is in contradiction with Ref. ^j] where the canonical version of SMM pre- 
dicted a independent of the size of the sources. The discrepancy between the cited canonical 
results and the present microcanonical results focuses the attention on the importance of 
using appropriate models for describing finite systems. 



H. Slope dependence on the emitted fragment size 

Taking into account that isoscaling has been argued on grandcanonical basis and above 
presented microcanonical results suggest important deviations from grandcanonical expec- 
tations, it would be interesting to check whether isoscaling is valid also for smaller sources. 

The results obtained for ((40, 20) and (50, 20)) at E ex = 5 MeV/nucleon and V = QV 
proved a linear behavior of the yield ratios with the isotopic composition of the emitted 
fragment together with a monotonic increase of the absolute value of the slope parameters 
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FIG. 5: Upper panels: Isotopic yield ratios for the primary decay of (50, 20) and (40, 20) at 
E ex = 5 MeV/nucleon and V = 6Vq as a function of N (left panel) and Z (right panel). Lower 
panels: Dependence of the slope parameters a and (3 of the size of the considered emitted cluster. 
Here a (/?) represents the slope of log R2%(Z, N) versus N (Z) for each considered value of Z (N). 
The considered sequence of N(Z) is 1,2,3,4,5,6,8,10. 

with the size of the emitted fragment as displayed in Fig. |5] where we plot the results 
obtained for fragments with 1 < Z < 10 and 1 < N < 10. This means that for small sources 
we obtain isoscaling in the sense of Eqs. (J7J) and (jSJ) without the constancy of the slope 
parameters with the size of the cluster. 

This result makes us go back to large fragmenting sources in order to evidence the de- 
pendence of a and (3 on the size of the emitted cluster. From our analysis results that by 
considering large emitted fragments isoscaling is obtained with slope parameters dependent 
on the size of the cluster. To illustrate this result, in Fig. Elwe plot cn(Z) for ((200, 80) and 
(185, 80)) at V — 6V and E ex = 5 MeV/nucleon. One can see that while for Z < 10 a is 
constant, for larger emitted fragments a increases monotonically with Z. It is interesting to 
mention that isoscaling slope parameters dependent on the emitted cluster size have been 
evidenced also in the framework of a dynamical stochastic mean field model [l^ . 

Because present detectors are not able to isotopically separate fragments with charge 
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FIG. 6: Variation of the slope parameter a with Z for the fragmenting systems (185,80) and 
(200,80) at V = GVq and E ex = 5 MeV/nucleon. Open symbols correspond to the primary stage, 
close symbols correspond to the asymptotic stage of the decay. 

larger than 8, the only case when one may hope to evidence this effect corresponds to light 
fragments emitted by light sources. 

The fragment dependence of the slope parameters can be interpreted as a finite size effect 
which, as expected, enhances with the diminish of the source and with the increase of the 
emitted fragment. 



IV. CONCLUSIONS 



To summarize, the present paper aims to investigate isoscaling in nuclear multifragmen- 
tation. In order to offer a complete and as accurate as possible image on the subject, the 
study was done using a microcanonical multifragmentation model and covered a large range 
of source sizes (A=40-200), excitation energies (E ex — 2-15 MeV/nucleon) and freeze-out 
volumes (V— 3Vo - 10Vb). Primary and asymptotic stages of the decay have been considered. 

The results indicate that despite in all considered cases the logarithm of the isotopic yields 
corresponding to the two different sources have a linear behavior as a function of Z and N, 
the slope parameters are strongly dependent on most observables characterizing the state 
of the sources. The only observable not affecting a and f3 is the freeze-out volume. This 
result could be of particular importance for experiments aiming to extract the asymmetry 
term of the binding energy using multifragmentation data, since most probably equilibrated 
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excited systems formed in heavy ion collisions have a distribution of volumes and the event 
separation as a function of volume is not by far possible. The isoscaling slope parameters 
are affected also by sequential particle emission. 

Deviations from strict isoscaling have been observed for heavy clusters originating from 
large multifragmenting systems and light clusters emitted by light systems and can be inter- 
preted as finite size effects. Such deviations from grandcanonical predictions are expected 
to manifest especially for small systems and stress the importance of using approximations 
suitable to the physical situation. 

Since the microcanonical ensemble is by principle the most appropriate statistical tool to 
describe the disintegration of small isolated systems and microcanonical multifragmentation 
models in general proved to be able to describe with accuracy experimental data, the most 
reliable description of isoscaling is expected to be obtained with this kind of models. The ac- 
curacy of isoscaling parameters is essential for obtaining correct information on fundamental 
quantities like the asymmetry term of the binding energy. 
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